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Wave  Momentum  Flux  Parameter 
for  Coastal  Structure  Design 

by  Steven  A.  Hughes 


PURPOSE:  The  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  described  herein 
provides  information  about  a  new  wave  parameter  for  characterizing  wave  processes  at  coastal 
structures.  A  description  of  the  parameter  is  given  along  with  sample  calculations  for  periodic  waves 
and  solitary  waves.  The  first  application  of  this  physically  relevant  parameter  has  been  development 
of  new  empirical  relationships  for  irregular  wave  runup  on  smooth,  impermeable  slopes  described  in 
CHETN-III-68  (Hughes  2003). 

COMMON  WAVE  PARAMETERS:  Coastal  engineers  have  established  useful  design  guidance 
by  augmenting  theoretical  reasoning  with  empirical  coefficients  determined  from  small-scale 
laboratory  testing.  Waves  are  usually  included  in  empirical  design  relationships  via  one  or  more 
wave  parameters  considered  to  be  representative  of  the  incident  wave  condition.  Common  regular 
and  irregular  wave  parameters  are  listed  in  Table  1. 


Table  1 

Common  Wave  and  Fluid  Parameters 

Regular  Wave  Parameters 

H  =  wave  height 

Hp  =  deepwater  wave  height 

L  =  local  wavelength 

Lp  =  deepwater  wavelength 

T  =  wave  period 

k  =  wave  number  [=  2n/L] 

Irregular  Wave  Parameters 

Hmo  =  zeroth-moment  wave  height 

Hs  =  significant  wave  height  [=  H1/3] 

=  root-mean-squared  wave  height 

Hw%  =  10%  of  waves  are  higher 

Tp  =  spectral  peak  wave  period 

7m  =  mean  wave  period 

Lp  =  wavelength  associated  with  Tp 

Lap  =  deepwater  wavelength  associated  with  Tp 

Lm  =  wavelength  associated  with  Tm 

Lom  =  deepwater  wavelength  associated  with  Tp, 

Fluid  and  Other  Parameters 

g  =  gravitational  acceleration 

p  =  fluid  density 

p  =  coefficient  of  dynamic  viscosity 

v  =  coefficient  of  kinematic  viscosity 

h  =  water  depth 

a  =  beach  or  structure  slope 

Sometimes,  these  wave  parameters  are  combined  to  form  dimensionless  variables  that  may  include 
relevant  fluid  and  other  parameters  such  as  those  given  in  Table  1 .  The  more  common  dimensionless 
wave  parameters  that  have  appeared  in  coastal  structure  design  guidance  are  shown  in  Table  2.  Note 
that  relative  wave  height  H/h  is  the  only  parameter  that  does  not  contain  wave  period  or  wavelength. 
Hence,  H/h  is  the  only  parameter  listed  that  is  applicable  to  solitary  waves. 

The  dimensionless  wave  parameters  listed  in  Table  2  strictly  pertain  to  uniform,  periodic  waves  of 
permanent  form  (i.e.,  regular  waves);  and  it  is  customary  to  use  first-order  wave  theory  to  calculate 
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values  for  wavelength.  These  same  dimensionless 
parameters  are  also  used  to  eharacterize  irregular 
wave  trains  by  substituting  wave  heights,  wave 
periods,  and  wavelengths  representative  of  irregu¬ 
lar  waves,  sueh  as  those  shown  on  Table  1. 

This  CHETN  describes  a  new  wave  parameter  for 
depicting  processes  that  occur  when  waves 
impinge  on  coastal  structures.  The  new  wave 
parameter  has  the  following  attributes: 

a.  The  parameter  is  physically  meaningful  so 
it  can  be  incorporated  into  theoretical 
models  of  specific  physical  processes  in  a 
rational  way. 

b.  The  parameter  applies  to  both  periodic  waves  and  nonperiodic  waves. 

c.  The  parameter  spans  the  range  of  relative  depths  from  deep  water  to  shallow  water. 

d.  The  parameter  is  stable  and  robust. 

e.  The  parameter  is  easy  to  estimate  so  design  guidance  using  the  parameter  can  be  pro¬ 
grammed  into  computer  spreadsheets  or  simple  programs. 

WAVE  MOMENTUM  FLUX  PARAMETER:  Surface  waves  possess  momentum  that  is  directed 
parallel  to  the  direction  of  wave  propagation.  When  a  wave  encounters  a  solid  object,  momentum  is 
reversed,  and  the  resulting  force  on  the  object  is  equal  to  the  rate  of  change  of  momentum,  also 
known  as  the  wave  momentum  flux.  Thus,  wave  momentum  flux  is  the  property  of  progressive 
waves  most  closely  related  to  force  loads  on  coastal  structures  or  any  other  solid  object  placed  in  the 
wave  field.  And  for  this  reason  wave  momentum  flux  is  a  compelling  wave  property  for 
characterizing  coastal  structure  response  to  wave  loading. 

The  instantaneous  flux  of  horizontal  momentum  {mj)  across  a  unit  area  of  a  vertical  plane  oriented 
parallel  to  the  wave  crests  is  given  by 

mf  (x,z,t)  =pd  +  pu  (1) 

where 

X  =  horizontal  direction  perpendicular  to  wave  crests 
z  =  vertical  direction,  positive  upward  with  z  =  0  at  still-water  level  (swl) 
t  ~  time 

Pd  =  instantaneous  wave  dynamic  pressure  at  a  specified  position 
p  =  water  density 

u  =  instantaneous  horizontal  water  velocity  at  the  same  specified  position 
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Integrating  the  instantaneous  wave  momentum  flux  over  the  water  depth,  i.e., 


Mp{x,t)=  ^(p^  +  pu^)dz  (2) 

J~h 

gives  the  total  wave  momentum  flux  per  unit  erest  length  acting  through  a  vertical  plane  parallel  to 
the  crest.  If  the  depth- integrated  Mp  is  also  integrated  over  the  wavelength,  the  result  becomes  the 
wave  “radiation  stress”  defined  by  Longuet-Higgins  and  Steward  (1964).  However,  there  is  sig¬ 
nificant  variation  of  depth-integrated  wave  momentum  flux  over  a  wavelength  from  large  positive 
values  to  large  negative  values.  So  instead  of  adopting  a  mean  value  which  is  quite  small  compared 
to  the  range  of  variation,  it  is  logical  when  considering  the  wave  force  loading  on  structures  to  focus 
on  the  maximum,  depth-averaged  wave  momentum  flux  that  occurs  during  passage  of  a  wave,  i.e., 
the  maximum  of  Equation  2,  which  occurs  at  the  wave  crest. 

Using  Equation  2,  maximum  depth-integrated  wave  momentum  flux  (MF)max  can  be  determined  for 
any  surface  wave  form  provided  the  velocity  and  pressure  field  can  be  specified.  Therefore,  {MF)max 
has  the  potential  of  being  a  unifying  wave  parameter  applicable  to  both  periodic  and  transient  wave 
types. 

ESTIMATING  MAXIMUM  DEPTH-INTEGRATED  WAVE  MOMENTUM  FLUX:  Hughes  (in 
preparation)  derived  formulas  for  estimating  the  maximum  depth-integrated  wave  momentum  flux 
for  periodic  (regular)  waves  and  solitary  waves. 

Periodic  Waves 


Expressions  derived  for  maximum  depth-integrated  wave  momentum  flux  from  linear  wave  theory 
do  not  include  that  part  of  the  wave  above  the  swl  where  a  significant  portion  of  the  wave  momen¬ 
tum  flux  is  found.  Estimates  of  (Mfymax  improve  using  extended-linear  theory  in  which  expressions 
for  linear  wave  kinematics  are  assumed  to  be  valid  in  the  crest  region.  However,  the  wave  form  is 
still  sinusoidal  rather  than  having  peaked  crests  and  shallow  troughs  typical  of  nonlinear  shoaled 
waves;  and  consequently,  the  theory  underpredicts  momentum  flux  under  the  crest.  Hughes  (in 
preparation)  calculated  values  of  (Mfymax  for  a  wide  range  of  nonlinear  uniform  waves  using  a 
numerical  technique  (Eourier  approximation)  that  optimized  the  solution  to  provide  the  best  fit  of  the 
fully  nonlinear  free  surface  boundary  conditions.  Results  were  expressed  in  terms  of  the 
dimensionless  parameter 


f 

V 


Mp 

pgh^ 


\ 

/max 


This  parameter  represents  the  nondimensional  maximum  depth- integrated  wave  momentum  flux,  and 
it  is  referred  to  as  the  “wave  momentum  flux  parameter.” 

Eigure  1  presents  the  dimensionless  wave  momentum  flux  parameter  versus  relative  depth  h/gl^. 
The  solid  lines  are  lines  of  constant  relative  wave  height  H/h.  Eor  a  constant  water  depth,  wave 
period  increases  toward  the  left  and  decreases  to  the  right.  The  range  of  relative  depths  covers  most 
coastal  applications. 
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Maximum  Wave  Momentum  Flux  -  Nonlinear  Theory 
Lines  of  Constant  H/h 


Figure  1 .  Dimensionless  wave  momentum  flux  versus  h/gf^ 


The  dashed  line  in  Figure  1  gives  the  steepness-limited  wave  breaking  criterion  tabulated  by 
Williams  (1985)  and  expressed  by  Sobey  (1998)  as  the  rational  approximation 


g 


=  Co  tanh 


V 


a^r 

\  +  byr  +  ^ 


(3) 


where 

Co  =  1.0575 

ai  =  0.7879 
2 

r  =  CO  h/g 
02  =  2.0064 
as  =  -0.0962 
bi  =  3.2924 
b2  =  -0.2645 

Sobey  noted  the  previous  expression  has  a  maximum  error  of  0.0014  over  the  range  of  Williams’ 
table.  Williams’  (1985)  tabulation  of  limit  waves  is  more  accurate  than  the  traditional  limit  wave 
steepness  given  by 
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—^  =  0.142  tanh(ytA)  (4) 

L 

which  overestimates  limiting  steepness  for  long  waves  and  underestimates  limiting  steepness  for 
short  waves. 


An  empirical  equation  for  estimating  the  wave  momentum  flux  parameter  for  finite  amplitude  steady 
waves  was  established  using  the  caleulated  curves  of  constant  H/h  shown  in  Figure  1 .  A  nonlinear 
best-fit  of  a  two  parameter  power  curve  was  performed  for  eaeh  calculated  ////i-curve,  and  the 
resulting  coefficients  and  exponents  for  each  power  curve  were  also  approximated  as  power  curves. 
The  resulting,  purely  emp ideal,  equation  representing  the  curves  of  constant ///A  shown  on  Figure  1 
is  given  as 


r 

V 


Mp 

pgh^ 


^  h  ^ 


-A, 


where 


do  =  0.639 


f— 1 

[h  J 


2.026 


4  =  0.180 


\n  j 


-0.391 


(5) 


(6) 

(V) 


Goodness-of-fit  of  Equation  5  compared  to  the  computed  values  given  in  Figure  1  is  shown  in 
Figure  2.  For  smaller  values  of  nondimensional  (Mf)max,  there  is  reasonable  correspondenee  except 
for  the  left-most  points  of  each  curve  (shown  below  the  line  of  equivalenee).  This  divergence  was 
caused  by  the  power  curve  tending  toward  positive  infinity  as  h/gT  0.  Greater  deviation  begins  to 
occur  for  dimensionless  (MF)max  >  0.6.  Nevertheless,  the  only  poorly  fitted  curve  is  for  H/h  =  0.8 
which  is  at  or  slightly  above  the  limiting  steepness  for  waves  on  a  horizontal  seabed.  This  poor 
correspondence  was  likely  the  result  of  foreing  the  numerical  computation  beyond  appropriate  limits. 

The  empirieal  equation  represented  by  Equation  5,  along  with  Equations  6  and  7,  provides  an  easy 
method  for  estimating  the  nondimensional  wave  momentum  flux  parameter  for  finite  amplitude, 
regular  waves.  However,  most  eoastal  structure  design  guidanee  developed  in  the  past  20-25  years 
has  used  wave  parameters  representative  of  unidirectional  irregular  wave  trains.  It  is  recommended 
that  the  wave  momentum  flux  parameter  for  irregular  wave  trains  be  estimated  by  substituting 
irregular  wave  parameters  Hmo  (zeroth-moment  wave  height)  and  Tp  (peak  spectral  period)  directly 
into  the  empirical  Equations  5,  6,  and  7.  Application  is  simple,  and  estimates  of  maximum  wave 
momentum  flux  should  be  reasonably  representative  of  the  irregular  wave  train.  However,  other 
irregular  wave  parameters  might  be  used  depending  on  the  applieation.  Therefore,  before  applying 
any  design  guidance  using  the  wave  momentum  flux  parameter,  it  is  important  to  ascertain  which 
wave  parameters  were  used  to  establish  that  particular  empirical  relationship. 
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Estimated  vs.  Calculated  Nonlinear  Mp 
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Figure  2.  Goodness-of-fit  of  nonlinear  momentum  flux  empirical  equation 
(Fourier  wave  theory) 


Solitary  Waves 

An  analytical  derivation  for  the  wave  momentum  flux  parameter  for  solitary  waves  was  given  by 
Hughes  (in  preparation)  as 


V 


M, 

pgh^ 


\n  j 


+ 


H 


\ 


+ 


Af2 

(H 

f 

r  M 

— +  1 

<  tan 

— +  1 

2M 

Vh  ; 

1 

L  2 

Vh  )_ 

1  3 

+  —  tan^ 
3 


MfH 

TtT 


(8) 


The  parameter  is  a  function  of  only  relative  wave  height  H/h  and  two  coefficients  M  and  N  that  are 
typically  presented  in  graphical  form.  To  accommodate  calculations,  a  nonlinear  curve  fit  was 
applied  to  the  plotted  M  and  A  curves  to  produce  the  following  simple  equations  that  give  reasonable 
values  for  M  and  N. 
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M  =  0.98 


2.24 


f— 1 

\h  J 


0.44 


(9) 


N  =  0.69  tanh 


2.38 

[.h  )_ 

(10) 


The  empirically-fit  equations  (solid  lines)  are  plotted  along  with  the  data  points  taken  from  the  Shore 
protection  manual  (1984)  in  Figure  3.  The  variation  of  maximum  nondimensional  wave  momentum 
flux  for  solitary  waves  as  a  function  of  H/h  is  shown  in  Figure  4.  These  values  represent  the  upper 
limit  of  nonlinear  (Fourier)  waves  when  h/{gT^  )  approaches  zero. 

APPLICATIONS:  It  is  anticipated  that  the  wave  momentum  flux  parameter  will  prove  useful  for 
developing  improved  semiempirical  formulas  to  describe  such  wave/structure  processes  as  runup, 
overtopping,  reflection,  transmission,  and  armor  stability.  As  an  example,  a  simple  expression  for 
wave  runup  on  plane,  impermeable  slopes  was  derived  by  assuming  the  weight  of  water  contained  in 
the  volume  above  the  still-water  line  at  maximum  runup  is  proportional  to  the  maximum  depth- 
integrated  wave  momentum  flux  in  the  wave  as  it  passes  the  structure  toe.  Hughes  (in  preparation) 
confirmed  this  simple  hypothesis  by  establishing  empirical  runup  relationships  based  on  the  wave 


Solitary  Wave  Theory 

Empirical  Fit  to  N  and  M  Functions 


Figure  3.  Nonlinear  curve  fit  to  solitary  wave  N  and  M  coefficients 
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Maximum  Wave  Momentum  Flux 

Solitary  Wave  Theory 


Figure  4.  Solitary  wave  nondimensional  wave  momentum  flux  versus  H/h 


momentum  flux  parameter  using  existing  runup  for  regular  waves  and  solitary  waves.  CHETN-III-68 
(Hughes  2003)  extended  this  concept  to  irregular  wave  runup,  providing  improved  estimation 
techniques  for  runup  on  smooth,  impermeable  slopes.  In  all  runup  cases,  structure  slope  was  found  to 
have  less  influence  than  previously  thought  for  the  range  of  slopes  in  the  data  set. 

EXAMPLE  CALCULATIONS:  The  following  examples  illustrate  calculation  of  the  wave 
momentum  flux  parameter  for  regular  and  solitary  waves. 


Example:  Regular  Wave  Momentum  Flux  Parameter 


Find:  The  nondimensional  wave  momentum  flux  parameter  and  the  corresponding  maximum  depth- 
integrated  wave  momentum  flux  per  unit  length  of  wave  crest  for  the  specified  regular  wave. 


Given: 

h  =  15  ft 

-  Water  depth 

7  =  8  s 

-  Regular  wave  period 

7/  =  9  ft 

-  Regular  wave  height 

g  =  32.2  ft/s^ 

-  Gravitational  acceleration 

pg  =  64  lb/ft3 

-  Specific  weight  of  sea  water 
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Estimate  Limit  Wave  Height:  First  check  to  be  certain  the  specified  wave  height  does  not  exceed 
the  steepness  limited  wave  height  for  the  given  wave  period  and  water  depth.  Calculate  the  circular 
wave  frequency  and  nondimensional  depth  as 


Ik 
(8  s) 


=  0.7854  1/s 


and 


r  = 


(0.7854  1/s)2(15  ft) 
(32.2  ft/s2) 


0.287 


The  limiting  wave  height  is  given  by  rearranging  Equation  3  and  replacing  the  constants  with  the 
given  numerical  values,  i.e., 


H 


limit 


Wj 


(1.0575)  tanh 


^  0.7879r  +  2.0064r2  -  0.0962r^  ^ 
l  +  3.2924r-0.2645r2 


Substituting  for  the  variables  O),  g,  and  r  yields 


H, 


limit 


32.2  ft/s^ 
(0.7854  1/s)^ 


(1.0575) 


tanh 


0.7879(0.287)  +  2.0064(0.287)^ -0.0962(0.287)^ 
1  +  3.2924  (0.287)  -  0.2645  (0.287)^ 


or 


(55.2  ft)  tanh  (0.2023)  =  11.0  ft 

Because  H  <  Humu,  the  specified  wave  condition  is  not  steepness  limited. 

Calculate  the  Wave  Momentum  Flux  Parameter:  First  calculate  values  of  relative  wave  height 
and  relative  depth  as 

//  9  ft  h  15  ft 

“  15  ft  “  ~  (32.2  ft/s2)(8  s)2  “ 

Next,  find  the  values  of  the  coefficient  Aq  and  from  Equations  6  and  7,  respectively  as 
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^0 

Ai 


0.639 

0.180 


f— 1 

[hj 

f-] 

[hj 


2.026 


-0.391 


=  0.639 


=  0.180  (0.6)-° 


=  0.2270 

=  0.2198 


Finally,  the  nondimensional  wave  momentum  flux  parameter  is  caleulated  from  Equation  5  as 


f 

V 


Mp 

pgh^ 


max 


“^0 


f 


0.2270  (0.0073)'°^'‘^*  =  0.669 


Calculate  the  Maximum  Depth-Integrated  Wave  Momentum  Flux:  The  dimensional  value  of 
maximum  depth-integrated  wave  momentum  flux  per  unit  length  along  the  wave  erest  is  obtained 
simply  as 


=  ip  g)  (h^) 


pgh^ 


(64  lb/ft3)  (15  ft)2  (0.669)  =  9,634  Ib/ft 


Example:  Solitary  Wave  Momentum  Flux  Parameter 


Find:  The  nondimensional  wave  momentum  flux  parameter  and  the  corresponding  maximum 
depth-integrated  wave  momentum  flux  per  unit  length  of  wave  crest  for  the  specified  solitary  wave. 

Given: 


h  = 

15  ft 

-  Water  depth 

H  = 

9  ft 

-  Solitary  wave  height 

g  = 

32.2  ft/s2 

-  Gravitational  acceleration 

Pg  = 

64  lb/ft3 

-  Specific  weight  of  sea  water 

Calculate  the  Wave  Momentum  Flux  Parameter:  This  solitary  wave  has  the  same  relative  wave 
height  as  the  previous  example,  i.e.. 


H 

h 


9  ft 
15  ft 


=  0.6 


Using  this  value  of  H/h,  determine  the  coefficients  M  and  N from  Equations  9  and  10,  respectively: 


M  =  0.98  <  tanh 


2.24 


0.44 


\n  j 


=  0.98  {tanh  [2.24  (0.6)]}0-44  =  0.923 
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and 


N  =  0.69  tanh 


2.38 


\n  j 


0.69  tanh  [2.38  (0.6)]  =  0.615 


A  quick  check  of  the  plots  in  Figure  3  with  a  value  of  H/h  =  0.6  eonfirms  the  ealeulations. 

Substituting  values  for  H/h,  M,  and  N  into  Equation  8  yields  the  nondimensional  wave  momentum 
flux  parameter  for  the  speeified  solitary  wave,  i.e., 


pgh^ 


(0.6f +2(0.6) 


,  (^(0.6.1)  ^  tan 
2(0.923)^  ’  ’ 


(0.923) 


(0.6  +  1) 


1  3 

+  — tan 
3 


(0.923) 


(0.6  +  1) 


j 


0.78  +  0.3278 


tan (0.7384) +  ^tan^  (0.7384) 


1.161 


The  plot  shown  in  Figure  4  agrees  with  the  previous  ealeulation  for  a  value  of  H/h  =  0.6. 

Caution:  If  you  perform  the  previous  ealeulation  on  a  hand-held  ealeulator,  you  may  need  to 
eonvert  the  arguments  of  the  tangent  funetion  from  radians  to  degrees  before  taking  the  tangent. 

Calculate  the  Maximum  Depth-Integrated  Wave  Momentum  Flux:  Just  as  for  the  regular  wave 
example,  the  dimensional  value  of  maximum  depth-integrated  wave  momentum  flux  per  unit  length 
along  the  wave  erest  is  ealeulated  as 


(MF)max  =  {pg)  (h^) 


Mp 

pgh^ 


(64  lb/ft3)  (15  ft)2  (1.161)  =  16.718  Ib/ft 


The  value  of  (MF)max  for  the  solitary  wave  is  1 .74  times  greater  than  the  regular  wave  having  the 
same  value  of  H/h.  In  faet,  the  solitary  wave  is  the  limit  of  the  regular  wave  as  wave  period  beeomes 
very  long. 


SUMMARY:  This  CHETN  has  deseribed  a  new  parameter  representing  the  maximum  depth- 
integrated  wave  momentum  flux  oeeurring  in  a  wave.  Beeause  wave  momentum  flux  has  units  of 
foree  per  unit  erest  width,  it  is  a  physieal  deseriptor  of  wave  forees  aeting  on  eoastal  struetures.  For 
periodie  waves  an  easily  applied  empirieal  expression  was  given  to  estimate  nondimensional 
maximum  depth-integrated  wave  momentum  flux  as  a  funetion  of  relative  wave  height  and  relative 
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water  depth.  A  formula  was  also  given  for  estimating  the  wave  momentum  flux  parameter  for 
solitary  waves.  Example  ealculations  are  provided  for  both  wave  types. 

The  wave  momentum  flux  parameter  is  useful  for  developing  improved  semiempirical  formulas  to 
deseribe  sueh  wave/strueture  proeesses  as  runup,  overtopping,  refleetion,  transmission,  and  armor 
stability.  The  first  applieation  of  the  parameter  has  been  revised  design  guidance  for  irregular  wave 
runup  on  smooth,  impermeable  plane  slopes  as  detailed  in  CHETN-III-68  (Hughes  2003). 

ADDITIONAL  INFORMATION:  This  CHETN  is  a  product  of  the  Scour  at  Inlet  Structures  Work 
Unit  of  the  Coastal  Inlets  Research  Program  (CIRP)  being  conducted  at  the  U.S.  Army  Engineer 
Research  and  Development  Center,  Coastal  and  Hydraulics  Laboratory.  Questions  about  this 
technical  note  can  be  addressed  to  Dr.  Steven  A.  Hughes  (Voice:  601-634-2026,  Fax:  601-634- 
3433,  email:  Steven.A.Hughes@erdc.usace.armv.mil).  For  information  about  the  Coastal  Inlets 
Research  Program  (CIRP),  please  contact  the  CIRP  Technical  Leader,  Dr.  Nicholas  C.  Kraus  at 
Nicholas.C.Kraus@erdc.usace.armv.mil.  Beneficial  reviews  were  provided  by  Mr.  Dennis  Markle 
and  Dr.  Jeff  Melby,  Coastal  and  Hydraulics  Laboratory;  and  Mr.  John  Ahrens,  retired  Coastal  and 
Hydraulics  Laboratory  and  Sea  Grant. 
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